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How an autonomous quantum Maxwell demon can harness correlated information 
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University of New Mexico, Albuquerque, NM 87131, USA 

We study an autonomous quantum system which exhibits refrigeration under an information-work 
trade-off like a Maxwell demon. The system becomes correlated as a single “demon” qubit interacts 
sequentially with memory qubits while in contact with two heat reservoirs of different temperatures. 
Using strong subadditivity of the von Neumann entropy, we derive a global Clausius inequality to 
show thermodynamic advantages from access to correlated information. It is demonstrated, in a 
matrix product density operator formalism, that our demon can simultaneously realize refrigeration 
against a thermal gradient and erasure of information from its memory, which is impossible without 
correlations. The phenomenon can be even enhanced by the presence of quantum coherence. 


I. INTRODUCTION 

With the discovery of his famous gedanken experiment, 
James Clerk Maxwell started physics on the long road to 
the unification of thermodynamics with information the¬ 
ory [1-3]. Nearly 200 years later, with the ever-advancing 
miniaturization of devices [4-8], the idea that one could 
exploit detailed knowledge of a system’s microstate as a 
thermodynamic resource no longer seems like a mere ide¬ 
alization [9-15]. To the contrary, we hope to someday 
soon push our technology to the regime where quantum 
effects become relevant [16-27]. Since there are more del¬ 
icate constraints on correlations in the quantum world, 
like strong subadditivity of von Neumann entropy, as well 
as a variety of monogamy relations (exclusive trade-offs) 
among different kinds of correlations, it is timely to con¬ 
sider a question unforeseen in Maxwell’s time: How can 
correlations be harnessed for the performance of these de¬ 
vices in the fully general quantum setting? 

This significant question has been studied extensively 
from different perspectives, and it has been mentioned 
that correlations can yield an advantage in performing 
certain thermodynamic tasks, such as work extraction 
[ .8, 28-35]. In this paper, we focus on two key chal¬ 
lenges: (i) accounting for the contribution of correlations 
to the second law of thermodynamics and (ii) comparing 
quantum and classical correlations as thermodynamic re¬ 
sources on the same footing. Regarding (i), it is impor¬ 
tant to notice that the system is coupled to external de¬ 
grees of freedom in quantum open-system dynamics, and 
the resulting correlations can be utilized, in principle, 
by an external, classical agent to “measure” the system 
and gain knowledge in the framework of feedback control 
[13]. However, this utilization of correlations should be 
distinguished from the way an internal agent like a quan¬ 
tum Maxwell demon, or a part of the quantum system, 
handles correlations, since the latter is the scenario we 
are interested in here. Regarding (ii), the majority of 
prior works [ .6, 19, 23, 24, 33, 36, 37] defines quantum 


and classical Maxwell demons in such a manner that the 
classical version is treated as a limiting case (typically 
with a fixed basis). This obviously makes the classical 
demon at most as powerful as the quantum one, but it 
may not immediately imply that classical correlations are 
at most as useful as quantum correlations for the demon 
who can handle both kinds. 

To address point (i), we study a quantum system which 
exhibits an information-work trade-off independently of 
any external agent. Following a pioneering work of the 
classical autonomous Maxwell demon studied in Ref. [9], 
we consider a scenario in which a qubit, in contact with 
two heat reservoirs of different temperatures, accesses se¬ 
quentially many memory qubits, and we analyze their 
correlations in the fully quantum regime. We think this 
to be among the simplest models to exhibit the ther¬ 
modynamic features akin to a quantum Maxwell demon. 
For point (ii), as defined formally in the next section, we 
stress that our demon is identical regardless of the na¬ 
ture of information in its memory. This allows us to un¬ 
ambiguously define work performed by the demon, since 
there has otherwise been difficulty in doing so due to the 
intrinsic “uncertainty” of quantum systems [38-47]. 

Here we discover a new nonequilibrium phase of our 
demon: refrigeration against a thermal gradient coinci¬ 
dent with memory erasure at the expense of correlations. 
We find that correlations enable our demon to exploit 
quantum coherence to realize an advantage over its clas¬ 
sical counterpart. To handle the complexity of correla¬ 
tions, we apply a tool from condensed-matter physics— 
the matrix product density operator formalism—which 
makes tractable our calculations in the fully quantum- 
correlated regime. This technique exploits the fact that 
correlations built under local interactions change locally 
[48-51]. It thus gives a powerful tool for treating correla¬ 
tions in the thermodynamic framework, a task which has 
otherwise remained formidable. We expect our proof-of- 
principle to inspire a cross-fertilization of quantum many- 
body physics and quantum information with quantum 
thermodynamics. 
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Section II is an introduction and overview of the au¬ 
tonomous quantum demon. In Sec. Ill A, we apply tech¬ 
niques from quantum information theory to derive its 
effective second law constraint. We give an analytic solu¬ 
tion using our matrix product density operator formalism 
in Sec. IIIB. In Sec. IIIC, we demonstrate the existence 
of a phase of simultaneous refrigeration and erasure in 
the presence of correlations, and in Sec. HID, we demon¬ 
strate the quantum advantage. Finally, we close with a 
discussion in Sec. IV. 


II. OVERVIEW OF THE MODEL 

Our autonomous Maxwell demon [Fig. 1(a)] is a gener¬ 
alization of the model proposed in Ref. [9] . The “demon” 
D is a qubit spanned by orthonormal basis states | g) and 
|e) with energy gap E e —E g = A. It interacts sequentially 
with qubits in an infinitely long one-dimensional array 
M, called the memory. The qubits in M are energeti¬ 
cally degenerate and may be initially correlated, though 
they are noninteracting among themselves. The memory 
thus acts only as an informational, and not an energy, 
resource. When discussing a given interaction, we will 
refer to the interacting memory qubit simply as M, the 
system of memory qubits which have previously inter¬ 
acted as M, and that of those which are yet to interact 
as M. For each such qubit, we define a “classical basis” 
{|0), 11)}. In what follows, we choose the qubit Pauli-z 
operator of D, cr z D , to be diagonal in its energy eigenbasis 
and that of M, cr z M to be diagonal in its classical basis. 
We are free to choose a phase convention, which we keep 
fixed, in defining the Pauli-x and -y operators, a x and 
a y respectively, of these systems. We will also need 

2d 0 = 1 + g z 2<j + = & x + ia y 
2<j 1 = 1 — a z 2<t _ = <j x — ia y 

on these systems as well. Finally, we will use 

C = (<5m) - (om> = (<*m) 

as a shorthand for the population bias of M in its z basis. 

The device operates cyclically; D interacts with M for 
an interaction time r before M moves by one site to the 
right, and the sequence repeats. Each interaction con¬ 
sists of two simultaneous processes, each in contact with 
a different thermal reservoir. In the first, the demon un¬ 
dergoes intrinsic transitions , wherein it exchanges a unit 
A of energy with a “hot” reservoir at temperature In 
the second, the joint system DM undergoes cooperative 
transitions , wherein the demon exchanges A of energy 
with a “cold” reservoir at temperature T c < X^, and the 
interacting qubit is flipped. Intrinsic transitions occur 
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FIG. 1. (Color online) (a) A snapshot of our quantum 
Maxwell demon. The demon qubit D interacts sequentially 
with each qubit in a memory M via an open-system process 
in contact with two thermal reservoirs of different tempera¬ 
tures. For a given interaction, M is the currently interacting 
qubit, M the subsystem of previously interacting qubits, and 
M that of qubits yet to interact. The joint state of DM. is 
described by a matrix product density operator with periodic 
boundary conditions, though we take sufficiently many inter¬ 
actions that the system has reached a periodic steady state. 
Interaction with the shared system D allows the possibility 
to build further correlations within M, hence the double line, 
(b) The transition diagram of an individual interaction be¬ 
tween the demon qubit and a memory qubit. Solid arrows 
represent possible stochastic transitions, allowed by exchang¬ 
ing a unit A of energy with the associated reservoir, with 
transition rates satisfying the detailed balance condition. For 
example, intrinsic transitions coupled with the thermal reser¬ 
voir of temperature Th change only the state of the demon, 
whereas cooperative transitions coupled with that of T c (< Th) 
change the joint state of the demon and the memory qubit, 
leaving a record on the memory. The dashed arrows indicate 
the cyclic processes which refrigerate (pump energy against 
the gradient), with the horizontal dot-dashed transitions pro¬ 
vided by the memory shift. Note both refrigerative processes 
result in a net flip of the memory qubit from “0” to “1.” 


with rates T g ^ e and T g< - e , and cooperative transitions 
with rates T^o-^ei and r^o^-ei* These rates are chosen to 
satisfy detailed balance 
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where the /% = 1/Ti are the inverse temperatures, and 
we have chosen units such that the Boltzmann factor is 
1. We thus describe each transition by a Lindblad jump 
operator 


— \/r -g—±e&]j ® Ijlf 
Lg^ e = yjT g ^ e d+ <S) 1m 


-LgO^-el — \/^ gO^-el^ d 
L gO-^el ^/^gOi—el^jj 


°M- 


Finally, we define 


e = tanh 


{Pc-Ph) A 

2 


as a parameter quantifying the magnitude of the thermal 
gradient. 

Formally, our interaction sequence is generated by the 
time-dependent Lindbladian 


Comit) = ^2 ^DM Q (j T ~ t)® [* ~ U ~ i) 7- ]. 

3 = 1 


where is the time-independent interaction Lindbla¬ 
dian between D and the jth qubit in M, described by 
the aforementioned Lindblad jump operators (see Ap¬ 
pendix A for details). 0 is the Heaviside step function, 
and so these factors “switch on” the interaction C^ M for 
t G [(j — 1)t, jr). In the limit of many interactions, the 
system reaches a periodic steady state for which 


Pdm = Pdm (nr) = p DM [(n + 1 )r], 


( 1 ) 


where 


Similarly, we will refer to the corresponding interaction 
Lindbladian simply as dropping the qubit label for 

convenience. 

There are two special cases of this model which are of 
interest to us. The first is the classical case, in which 
the evolution during each interaction is constrained to 
population dynamics in the eigenbasis of cr z D (8) cc z M . The 
interaction Cdm is “classical” in the sense that if the 
initial state of DM is diagonal in this basis, then the 
evolution will remain so throughout. That is to say that 
the evolution does not mix together populations and co¬ 
herences in this basis. Additionally, it drives the state of 
DM to a unique fixed point p = p^ (8) p ^, which 
is a product of classical states of D and M. 

The second special case is the uncorrelated case, in 
which the state of M is initially a product, and we neglect 
any correlations that are built over repeated interactions 
with the demon. Though the full dynamics, generated 
by is a sequence of local quantum operations, it 

may either build or consume correlations in the full state, 
since there is a common degree of freedom, D, between 
all of them. In general, one or both of these special cases 
may hold. 

Figure 1(b) shows two cyclic processes over which en¬ 
ergy is pumped against the thermal gradient (dashed ar¬ 
rows). The horizontal bit-flip transitions (dot-dashed ar¬ 
rows), which complete each cycle, are provided by the 
translation of the memory. If the sum of the probabili¬ 
ties of these processes is greater than that for the reverse 
processes, then there is a net flow of energy against the 
gradient, and we say the system is refrigerating. Note 
that both processes flip the qubit from the |0) state to 
the |1) state, so a record of the refrigeration is left on the 
memory. We thus define 


Qh^-c 


A 

2~ 


(C'-C) 


as the heat flow from the hot reservoir to the cold reser¬ 
voir during the interaction, which is negative when the 
system is refrigerating. Note that this property of the 
model removes any ambiguity in our definition of work. 
We also define 


pDmif) = [Adm(O)] = e^o c ™(s)ds [p DM (0)] (2) 


A S m = S(p' m )-S(pm). 


is the state of the full system at time £, and n G Z + is a 
sufficiently large positive integer. In what follows, we will 
only be interested in the performance of the device over a 
single interaction in periodic steady state. We therefore 
omit the explicit time dependence where appropriate and 
denote quantities corresponding to the outgoing qubit in 
the interaction as primed and those to the incoming as 
unprimed {e.g., p M = Pm (jit ) and p’ M = p M [(n + l)r]}. 


where 

S(p ) = — tr (p log p) 


is the von Neumann entropy. The change A Sm rep¬ 
resents the entropy dumped onto the memory. When 
A Sm < 0, the memory is being erased. 
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In the uncorrelated classical regime, our model reduces 
to that of Ref. [ ], where it was first introduced. There, 
it was shown that 

Qh^c(Pc — Ph) + ASm > 0, (3) 

in this regime. We refer to Eq. (3) as the local Clau¬ 
sius inequality. It represents a strict trade-off between 
the refrigeration (Qh-> c < 0) and erasure (A Sm < 0) 
phases of the device. This trade-off can be understood 
geometrically. In the absence of correlations or quantum 
coherence, the interacting qubit is driven monotonically 
along the z axis of the Bloch sphere to its fixed point, 
which has positive bias = e > 0. Erasure and 
refrigeration correspond to approaching this fixed point 
from C < £ and £ > £, respectively. So long as this 
evolution is monotonic, the trade-off is necessarily strict. 

It is noted in Ref. [9] that when correlations cannot be 
neglected, this result generalizes to 

Qh^c(Pc ~ Ph) + A Sm > A Id-.M: (4) 

where 

A I A:B = I (A : By - I (A : B) 
is the change in the quantum mutual information 

I(A : B) = S(p A ) + S(p B ) ~ S(p AB ). 

A Id.m can be either positive or negative in general, so 
the trade-off in Eq. (3) becomes no longer strict. This 
suggests the possibility of a third thermodynamically 
nontrivial phase, in which the device is refrigerating and 
erasing simultaneously. We expect this to occur when 
correlations between D and M induce nonmonotonic evo¬ 
lution on M alone. However, Eq. (4) is inconvenient in 
that it does not involve only terms that can be calcu¬ 
lated from the evolution on the memory alone and so 
requires knowledge of the periodic steady state, Eq. (1). 
One of our contributions is to provide a constraint free 
of this limitation, which is also strict in the presence of 
correlations and makes their utility as a thermodynamic 
resource manifest. 


III. MAIN RESULTS 

A. Global Clausius inequality 

Our first result is to recover the second law in terms of 
the relevant thermodynamic variables in the fully quan¬ 


tum correlated case. This is the global Clausius inequal¬ 
ity , which we state as a theorem. 

Theorem 1. The global Clausius inequality, 

Qh^dPc Ph) + A Sm ~ A/ M:i g - (5) 

which represents an information-work trade-off of the 
thermodynamic quantities defined above, holds for a peri¬ 
odic steady state in the fully quantum correlated regime. 
Accordingly, simultaneous refrigeration and erasure is 
possible only if it is attended by the consumption of cor¬ 
relations. 

Equation (5) represents a strict three-way trade-off 
among refrigeration, erasure, and the generation of cor¬ 
relations A I M .y L - 1 all of which may be calculated from 
knowledge of the initial and final states of the memory 
alone. The minus sign on A I m .m places consumption of 
correlations as a resource; the more negative this term, 
the more negative the local terms can be. 

To provide a sketch of the proof of Theorem 1 (full 
details can be found in Appendix B), we first note 

Qh^-c(Pc ~ Ph) + (^) 

which follows in the same manner as Eq. (4) from mono¬ 
tonic evolution to the fixed point , but including the 
system M as an ancilla. Next, we have 

A Wm = I(D : MM)' - I(D : MM) 

> I(D : M)' - I(D : MM) 

= I(D : MM) - I(D : MM) = 0. 

The first inequality follows from strong subaddivity: Mu¬ 
tual information is nonincreasing under partial trace. 
The last step follows from the periodic steady-state con¬ 
dition, Eq. (1). That is, the mutual information between 
D and the qubits yet to interact— MM before the inter¬ 
action and M after the interaction—is a constant across 
interactions in steady state. Thus 

Qh^c{Pc ~ Ph) T AS M m > 0, 

and using 

: — A I M: mi 


where A= 0 due to the fact that M does not partic¬ 
ipate in the interaction, therefore gives the result. 
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We see that simultaneous refrigeration and local era¬ 
sure is not forbidden by Eq. (5) so long as correlations 
are consumed, but it remains to be seen that our model 
actually exhibits such behavior. We address this in the 
following sections. 


B. Matrix product density operator solution 

The matrix product state formalism has been very suc¬ 
cessful as an efficient representation of the quantum cor¬ 
relations present in a one-dimensional spin-chain [52-5 z ]. 
Because dynamics in contact with a thermal bath are 
dissipative, we need a generalization of this formalism to 
mixed states and classical correlations. This is the ma¬ 
trix product density operator (MPDO) formalism, intro¬ 
duced in Refs. [ 18, 49] and further developed in Refs. [55- 
57]. The MPDO description of a mixed state of N d- 
dimensional spins with periodic boundary conditions is 
given by 


( N \ N 

V A lj J (7) 

o =i J i=i 

where the ij G {0,1,+,—}, and the A 1 * are Xj-i x Xj 
matrices corresponding to the spin at site j. V is the 
path-ordering operator, which places operators of higher 
j to the left in the product. It was shown in Ref. [48] 
that the description Eq. (7) can be obtained from a cor¬ 
responding pure matrix product state by tracing out an¬ 
cillary degrees of freedom on each of the spins. It thus 
reduces to those for pure matrix product states and clas¬ 
sically correlated distributions as special cases. Because 
we are interested in the periodic steady-state behavior of 
our model, we restrict ourselves to the single-site trans- 
lationally invariant case where the A 1 * are independent 
of the site label j. We thus have that the Xj = X are 
all equal. This quantity is known as the bond dimen¬ 
sion of the MPDO, which quantifies the degree to which 
the state is correlated. Note that, in this representa¬ 
tion, the state is specified by d 2 x 2 complex parameters, 
an exponential improvement — when bond dimension is 
polynomial in TV — over the 0(d N ) complex parameters 
required to specify each of the matrix elements of Pmpdo 
individually. 

In our model, we let Pm(0) = Pmpdo be the 

parametrization of the initial quantum correlated state 
of the memory M, which is uncorrelated with the initial 
state of D. We first solve the Lindblad master equation 
for an individual interaction 


dp dm 
dt 


£dM ( PDM ) 


analytically in MATHEMATICA, obtaining an expres¬ 
sion for the quantum operation 



Pdm(0) ^ YLT 


PdmO 11 ) 


~(ss) 

Pdm 



FIG. 2. (Color online) A pictorial representation of our time 
evolution method, which relies on a matrix product density 
operator (MPDO) description, for n — 3. (a) On the left is 
a schematic depiction of our interaction sequence, with quan¬ 
tum operations vectorized to let us represent the sequence as 
a quantum circuit. The shared degree of freedom is the demon 
qubit D, and the remaining degrees of freedom are the mem¬ 
ory M. We indicate that a particular qubit is not acted upon 
by a given quantum operation by drawing the qubit’s line over 
the quantum operation’s box. We parametrize the sequence 
as a matrix product operator (MPO), shown on the right, 
with physical indices labeled according to the corresponding 
qubit degrees of freedom, (b) We multiply our initial MPDO, 
for which the state of M is translationally invariant and un¬ 
correlated with that of D, by the MPO describing our interac¬ 
tion sequence to obtain an approximation to the steady state. 
This is done by performing a sequence of index contractions, 
shown on the left as dotted boxes. The approximation to the 
steady state, which is valid for sufficiently large n, is shown 
as the resulting MPDO on the right. All boundary conditions 
are taken to be periodic. 


cf) T = e CDMT . 


We then update our initial state Pdim( 0) to pDm( nT ) ~ 
p^i according to the sequential evolution, Eq. (2), using 
the method schematically depicted in Fig. 2. Here we 
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have vectorized the Hilbert space, representing quantum 
operations as matrices and density matrices as vectors 
(see Ref. [49] for details). In the figure, we represent 
each on-site tensor A lj by a shape (box or circle) and that 
tensor’s indices (either ij or the matrix indices of A lj for 
a particular value of ij) as lines from the corresponding 
shape. Tensors with a contracted index between them are 
joined by a line, and uncontracted indices are represented 
by open lines. This pictorial scheme is thus a convenient 
representation for tensor contractions. 

In Fig. 2(a), we apply a generalization of the re¬ 
sult from Ref. [58] to parametrize our series of sequen¬ 
tial quantum operations as a matrix product operator 
(MPO). Figure 2(b) shows the update, which is per¬ 
formed by simply multiplying the initial MPDO, Pdm{ 0), 
by the MPO, <F nr . This scheme permits us to only store 
3 x 4y 2 parameters: the tensor corresponding to MM, 
that to M, and a boundary tensor to D. The full details 
are given in Appendix C. 

We next apply the classicality of our interaction, </> T , 
to simplify our MPDO description for t erms 

of operationally defined quantities as our second main 
result. 

Theorem 2. The periodic steady state of DMM prior 
to the interaction on M may be expressed as 


A I m -m> same values for the input states pm{ 0) 

and [pm(0)] in steady state, where 

^(m)=(e-^ 

constitutes an individual rotation about the z axis applied 
transversally to every qubit in M. 

To see this, we first note that we calculate ~ for 
[/4m (0)] by making the replacement 

AT) (V 7/*,<*> (\ 

' MM MM \ MM) 

in Eq. (8). This is because transversal phase rotations 
affect neither the probabilities of classical measurements 

Pj: nor the states p$ by their definition. Furthermore, 

the states p^ are classical and so invariant under phase 
rotation. We thus have 

Y,PkP { D ®U Z mm ( Pmm ) = Ep* Kmm ( P ( d®Pmm ) 

k k 

= ^DMM ( PDMM ) • 


P { DMM = Y.PkP { D ®Pmm’ (4 

k 

where the k G {0, l} n are classical bit strings of length 
n, which correspond to classical records on the mem¬ 
ory subsystem M. p^ is the probability of the record k 

C^kPk = ^)> an d P^mm ™ hh e reduced state of MM con- 

'-(k) 

ditioned upon that record. pr D is the reduced state of D 
upon interacting with the string of uncorrelated classical 
pure qubits specified by k, which is classical. 

The proof and formal expressions for p £, p^\ and 
can be found in Appendix C. Theorem 2 says that the 
steady state is always separable and given by a convex 
combination over all possible classical histories of the 

memory. Note that the memory part P M] ^ could contain 
quantum coherence. It is convenient in that it expresses 
the steady state P^ M] ^ i n terms of the input (the p^ and 

the P^m) an d interaction (the p$) separately. This 
will allow us to engineer states of M to achieve a desired 
thermodynamic behavior. It is also a useful analytic tool 
when it is tractable (i.e., there are relatively few terms 
in the sum). We further demonstrate its utility by using 
it to prove that the demon’s performance in steady state 
is invariant under local phase rotations on the memory, 
which we state formally in the following corollary: 

Corollary 2.1. The Clausius terms, Qh^-c, ASm, and 


That is 

(I D ® U z "){p% s) M ) = WdUPdA (9) 

where 1 is the identity superoperator. We see from 
Eq. (9) that Cdm commutes with on Pp%, as the 
cooperative transition term commutes with and 

the intrinsic transition and Hamiltonian terms commute 
with Um? (Appendix A). Corollary 2.1 therefore follows 
from the fact that all of the terms in Eq. (5) are invariant 
under transversal phase rotations on p M1 ^ and p' MI ^- 

C. Simultaneous refrigeration and erasure 

We first examine a special case, which is simple enough 
that the calculation of Eq. (8) is analytically tractable 
but also correlated enough to reveal nontrivial thermo¬ 
dynamic behavior. 

Observation 1. Consider the family of Greenberger- 
Horne-Zeilinger- (GHZ) correlated states parametrized by 
a bias (, 

Mm = (yrrcio) 0iv +ywii) 0iv ), (io) 

with N sufficiently large to allow the system to reach a pe¬ 
riodic steady state. There exists a nonequilibrium phase 
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of simultaneous refrigeration (Qh^c < 0/ and erasure of 
memory (ASm < 0), when the memory is initially pre¬ 
pared to be the quantum entangled state \ip) M with a bias 
( and a temperature gradient e in the region illustrated 
in Fig. 3(a). 

We see that the reduced state on each qubit for any state 
in the family Eq. (10) is locally classical with bias £ and 
that for ( = =bl, IV^m is a product. For £ G (—1,1), how¬ 
ever, the state is entangled, with maximal entanglement 
at £ = 0. 

For this family, we calculate the Clausius terms in 
Eq. (5) analytically by letting the system undergo n in¬ 
teractions — so it reaches periodic steady state — and 
tracing out M. In Eq. (8), only two classical histories 
appear in the sum, and we have 


Pmm ~ / > 

/cG{0,l} L 


i + (-i) fe C 




( 11 ) 


and 


Pmm 

/ce{o,i} L 


l + (-i) fc C 


Pm ® (**) 




( 12 ) 


where 


PM ^ — tr DM |^ > (n+l)r Pd{ 0) 0 &M 



is the state of M following its interaction for an uncorre¬ 
lated tape of pure qubits in state |fc), in periodic steady 
state. The global entropy change AS M] ^ then takes the 
particularly simple form 


^MM ~ 2 


(1 


C) S(p'^) + (1 - C)S(p'^) 


Note that this quantity is always non-negative. By com¬ 
parison, the local quantities 


AS m = S 



-1(0) 

Pm 


-S 


^ (i + (& z ) 



and 


A 

~2 



i-c 

2 


yi) 


-C 





(a) 



FIG. 3. (Color online) (a) The nonequilibrium phase diagram 
for input states from the GHZ-correlated family described in 
Observation 1, analogous to Fig. 2 in Ref. [9], for r = 0.3 
(units of 1/A for h— 1). For f = ±1, the input is a product, 
and it is entangled for £ G (—1,1), with the strongest corre¬ 
lations at £ = 0, for which it is the GHZ state. Note that, 
for every state in this family, the reduced state of any qubit 
is diagonal in the z basis. We see that correlations shift the 
phase boundaries to induce a region of simultaneous refrigera¬ 
tion and erasure near the uncorrelated phase transition triple 
point (e = C = 0). (b) Dependence of the Clausius terms, 
Qh^c(/3 C — Ph) and ASm , on the interaction time, r when 
e — 0.01 in comparing the GHZ state for £ = 0—solid black 
and orange (gray) curves, respectively—to the product of its 
single-qubit reduced states, which are all maximally mixed— 
dotted black and orange (gray) curves, respectively. The re¬ 
duced evolution of the interacting memory qubit is monotonic 
when there is a strict trade-off of refrigeration and erasure; 
both cannot be simultaneously negative. This is obeyed by 
the uncorrelated input, but the correlated input shows a de¬ 
parture from monotonicity for r < 2.2. For long-enough in¬ 
teraction times, the effect of correlations is “washed out” as 
the terms for correlated input approach the values for their 
uncorrelated-input counterparts. 


Qh^-c 






































where = tr may be either positive or neg¬ 

ative. Using the expressions for the given in Ref. [9], 
we analytically construct an analogous nonequilibrium 
phase diagram, Fig. 3(a), for r = 0.3 (in units of 1/A, 
where we have set h = 1). We see that the effect of the 
correlation is to introduce a region of simultaneous refrig¬ 
eration and erasure near the uncorrelated phase transi¬ 
tion triple point at 6 = ( = 0. ( thus decreases below zero 
over the interaction in this region, indicating a deviation 
from monotonicity in the reduced evolution of M. This 
is the primary role of correlations in our thermodynamic 
model: in periodic steady state, initial correlations in M 
induce classical correlations between D and M prior to 
their interaction, allowing p M to pass through hitherto 
inaccessible regions of the Hilbert space in its evolution. 

We examine this more closely in Fig. 3(b), where 
we plot the quantities in the local Clausius inequality, 
Eq. (3), for the input given in Eq. (10) as a function of r 
for 5 = 0.01 and £ = 0 as compared to the corresponding 
terms for product input with the same reduced state on 
every qubit. In Fig. 4, we plot the generation of corre¬ 
lations A I m .m = A Sm ~ AS Mi g. We first note that, 
though correlations are always being consumed, the de¬ 
mon is not always able to use this to its advantage. This 
is as we expect, since the behavior must approach that of 
the uncorrelated case as r oo. For r < 0.54, we see in¬ 
creasing simultaneous refrigeration and erasure with ad¬ 
ditional interaction time. Past r ~ 0.54, additional inter¬ 
action time does not afford better thermodynamic perfor¬ 
mance, and, eventually, we see that correlations actually 
begin to hinder the demon’s erasure at r « 1.6. Finally, 
we see that the demon’s behavior approaches that of the 
uncorrelated case as r becomes large, as we expect. 

We see from the diagonal form of p MI[ that we would 
have had the same effect if the input had been prepared 
in the equivalent classically correlated (i.e., diagonal in 
the computational basis) state 



FIG. 4. (Color online) Generation of correlations A I M .p A be¬ 
tween the interacting qubit M and the qubits M yet to inter¬ 
act, as a function of the interaction time r, for the GHZ input 
of Fig. 3(b). r is in units of 1/A for h = 1, and £ = 0.01. Cor¬ 
relations are always being consumed, and so this term offsets 
the others in the Clausius inequality such that simultaneous 
refrigeration and erasure does not violate the second law for 
this input family. The dotted line is at A/ M:i q = — In 2 and 
is meant to provide a visual aid, as it represents the maximal 
possible consumption of correlations for this state. 

with 0 G [0,7r] ; and <fi G [0,27r), we define two states: 
(i) the GHZ-correlated input family in an ft basis p^ = 

IV’XV’Im su °h ^at 

I^)m = (\A + Cn I +fi)® N + \A - Cn I, 

(13) 

with N sufficiently large to allow the system to reach 
steady state, and (ii) its corresponding classically- 
correlated family, 


dc) _ 
Pm — 


1 - 
2 . 


(1 + C) (+“)*"+ (!-<) (a 1 ) 


] \ 0Af 


We are then led to ask: Is there anything to be said for 
the role of quantumness in this model? To address this, 
we consider another example, which will give our final 
result. 


D. Quantum thermodynamic advantage 

As a final result, we observe that correlations can en¬ 
able our model to exploit quantum coherence to gain a 
thermodynamic advantage 

Observation 2. Using an arbitrary orthonormal n basis, 
|+n) = cos (0/2) |0) + e l< t sin (0/2) |1), 

| —ft) = sin (0/2) |0) — e l(t) cos (0/2) |1), 


Pm = E ^0**'j Pm ( 14 ) 

with the kj G {0,1}, obtained by eliminating all off- 
diagonal matrix elements from pj^ in the z basis. The 
quantum entangled state p^ is advantageous in memory 

erasure over the classically correlated mixed state pj^ , 
whose bias for every single-qubit reduced state is identi¬ 
cal to that of p^p. 

This comparison might correspond to the scenario where 
the experimenter only has the ability to prepare correla¬ 
tions in a fixed basis, versus that where the experimenter 
has this ability in addition to the ability to perform phase 
rotations in an orthogonal basis. Thanks to Corollary 
2.1, we know that this second ability is the only addition 
over the first needed to observe the full range of thermo¬ 
dynamic performance. Because the states p^ and p^ 
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FIG. 5. (Color online) (a) The three-dimensional surface 
whose enclosure is the set of local Bloch vectors for which 
the locally coherent GHZ-correlated family attains a quantum 
erasure advantage over its corresponding classically correlated 
family, whose off-diagonal elements in the z basis have been 
projected out. Here r = 0.3 and £ = 0.01. We see clearly 
the predicted azimuthal symmetry about the z axis, and the 
advantage vanishes along this axis as we expect, since the two 
inputs are the same in this case, (b) The difference between 
the erasure performance for the quantum correlated locally 
coherent input state and that of its decohered counterpart for 
e — 0.01, 4> — 7t/2, 6 — 0.01, and Cn — —0.02 as a function of 
the interaction time r. A quantum advantage is demonstrated 
when this quantity is negative. We see that the advantage is 
several orders of magnitude smaller and lost much faster than 
the advantage afforded by classical correlations of Fig. 3(b). 


allow for all possible classical histories with some prob¬ 
ability, the demon’s performance cannot be easily calcu¬ 
lated with Eq. (8) as with Eq. (10). We thus obtain our 
results numerically using the MPDO description. 

In Fig. 5(a), we plot the exact region for which the ad¬ 
vantage exists [i.e., A< Aand A< 0]. We 
see here clearly the azimuthal symmetry of Corollary 2.1 


and that this advantage sharply disappears as 0 0. In 

Fig. 5(b) for e = 0.01, we show the difference in erasures 
between the quantum and classical inputs as a function 
of r. The quantum coherent state has an advantage 
when this quantity is negative. This plot bears some re¬ 
semblance to Fig. 3(b). This is no coincidence. In the 
same way as classical correlations in DM can allow the 
Bloch vector of M to momentarily move away from its 
fixed point along the z axis, they can also allow it to do 
so in the x-y plane, though the effect is smaller by several 
orders of magnitude. This deviation from monotonicity 
appears as an erasure advantage since our interaction </> T 
treats the classical and coherent components of the evo¬ 
lution separately. The refrigerative performances of the 
two inputs are thus identical as we expect. 

What is remarkable is that even though both pjj^ and 
p{j^ share the same classical histories in Eq. (8) , the small 

amount of coherence on the reduced state p^ seems to 
be able to seed the production of more coherence. This 
is especially surprising given the fact that our dynamics 
ct> r also has many classical features. The explanation is 
that the coherences on the conditional states of M for 
each classical history interfere such that there is more 
coherence on pM after the evolution than before. 


IV. CONCLUSION 

We have shown that our model of an autonomous 
Maxwell demon can utilize correlations in its memory 
to gain a thermodynamic advantage over its uncorre¬ 
lated counterpart. That is, as predicted by the global 
Clausius inequality, it performs simultaneous refrigera¬ 
tion and erasure, which is impossible without correla¬ 
tions, and exploits quantum coherence to enhance the 
performance of memory erasure. This work represents a 
proof-of-principle, and we hope that it paves the way to 
invent a quantum thermodynamic device which is more 
well suited to utilizing correlations. 

There are two concluding discussions. First, in this 
work, entanglement did not demonstrate a genuine ad¬ 
vantage over classical correlations. This is because of the 
long-range nature of GHZ correlations. Namely, trac¬ 
ing out a distant subsystem gives a classically corre¬ 
lated state but does not affect the thermodynamic per¬ 
formance, which is determined by the local interaction 
between the demon qubit and one of memory qubits. 
However, as far as we have explored, some translation- 
ally invariant, short-range-correlated states (such as the 
one-dimensional cluster state [59]) did not exhibit the 
simultaneous refrigeration-erasure phase. More demand¬ 
ing numerical calculation seems to be needed to gener¬ 
alize beyond these simple initial states, to address the 
question about strict supremacy of entanglement. 

Second, the simplicity of our model could bring the 
testing of quantum coherence as a thermodynamic re¬ 
source closer to experimental feasibility, since the model 
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shares similarities with quantum heat engines [60-7' ]. 
Though we consider a concrete experimental proposal to 
be beyond the current scope of the paper, we note that 
there are similarities between our model and the quan¬ 
tum heat engines considered in Refs. [72-74]. Here, a 
two-level system (TLS), such as a pair of electronic en¬ 
ergy levels of an optically active rubidium atom, is peri¬ 
odically driven by a red-detuned laser field and indepen¬ 
dently coupled to two thermal baths, such as the photon 
bath of a cavity and the “dephasing bath” of a buffer gas, 
at different temperatures. We envision the TLS serving 
the role of the demon in these systems, while the laser re¬ 
alizes the memory. As a buffer gas molecule approaches 
the rubidium atom, it perturbs the energy gap of the 
TLS into resonance with the laser, allowing the atom 
to become excited. Once the buffer molecule has scat¬ 
tered, and the gap of the TLS has widened back to its 
unperturbed value, the rubidium atom may fluoresce, ra¬ 
diating a quantum of energy into the photon bath. The 
authors of Ref. [74] note the existence of a regime where 
the dephasing bath may be cooled at the expense of work 
provided by the laser, with energy transfer to the pho¬ 
ton bath mediated by the TLS. This regime has been 
observed experimentally as reported in Refs. [7 2, 73]. 
Though we do not expect our model’s Hamiltionian to 
be isomorphic to that of this system, we do expect that 
it can be treated using our formalism in the case where 
the laser field state is correlated. 

The most severe imperfections to this setup would 
come from limitations on the experimenter’s ability to 
produce a precisely correlated state. We refrain from 
speculation on what experimental techniques are needed 
to prepare such states. Here we found that the novel ther¬ 
modynamic behavior was rather sensitive to the strength 
of the correlations [in the phase diagram, Fig. 3(a), the 
simultaneous-refrigeration-erasure phase occupies a nar¬ 
row range of £], but not so sensitive to classical corre¬ 
lations versus entanglement. This latter feature, should 
it survive in other architectures, could be promising for 
prospective experiments. Finally, our investigation was 
by no means exhaustive, so it remains to be seen what 
types of correlation are needed to optimize the thermody¬ 
namic performance of these particular experimental se¬ 
tups. It would be interesting to see if the utility of cor¬ 
relations as a thermodynamic resource survives in these 
architectures as well. 


Appendix A: Exposition of the Interaction 
Lindbladian 

Here we give the details of our interaction Lindbladian 
jCdm • As mentioned in the main text, the intrinsic and 
cooperative transitions are characterized by the Lindblad 
jump operators 


Ag—>. e D 

Lg-^e — D ® -fl-M 


for the intrinsic transitions, and 


LgO^-el — \/r g 0 _>-el (7 jj ® ® M 
LgO-^el el^D ® O’ M 

for the cooperative transitions. Here, <7 ± = \ ( a x ± ia y ) 
are the qubit ladder operators, and a^ x,y ^ are qubit Pauli- 
x and -y operators. Once again, the transition rates are 
chosen to satisfy detailed balance, 

_ e -p h A Fgp^el = e -/3 c A. 

IVe IWel 


Finally, the demon has an intrinsic Hamiltonian, 
Hd= A|e) (e\ = ^(i D -*b), 


and so the full dynamics of the interaction are generated 
by the Lindbladian 


&DM ( PDM ) = — i[Hd^M, PDm]+^ (PDM ) , (Al) 

i 


where 


( PDM ) = Lip DM L[ - -{L\Li,f)DM }j 
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is the term corresponding to a particular jump operator. 
[.,.] and {.,.} are the commutator and anticommutator, 
respectively, and we have chosen units such that h = 
1. As mentioned in the main text, C g o-^ e i and C g o^ e i 
commute with as defined in Corollary 2.1, due to 

the fact that a z D 0 o z M commutes with a l D 0 cr J M for i,j G 
{x,y}. Clearly C g ^. e , C g ^ e , and the Hamiltonian terms 
commute with . 
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jCdm has a fixed point p^p for which 



and given by a product 


from the periodic steady-state condition on D , together 
with the fact that the interaction acts as the identity on 
M, Eq. (Bl) reduces to 


AS 


DMM 


- tr 


(Pm ~ Pm) ^P { m 


> 0 


(B2) 


Next, we note that 




where 


s(/p) _ 1 ( 1 0 \ 

p D i + e p h A ^ 0 e PhA ) 


in the energy eigenbasis, and 

zifp) = _I_ ( 1 0 ^ 

Pm 1 + e (( 3 h -l 3 c )A ^ 0 e (^-/ 3 c)A y 


ln PM } ) = - &) (l-C)-ln 1 + e^ ^ c)A 


All of the constant terms cancel in the difference, and so 
we are left with 


tr 


(p'm - pm) In pti p) ] = |(c' - cm - Ph) 


Defining Q h ^ c = f (£' - C), Eq. (B2) gives 


Qh—>c(Pc Ph) + AS DMM — ® (E3) 


in the z basis. The form of this fixed point will be im¬ 
portant for the following proofs. 


Appendix B: Generalized Clausius Inequality 

Here we give a proof of Eqs. (3), (4), and (6). As 
stated in the main text, we have monotonic evolution to 
the fixed point in a single interaction (</> T (g)I^) 


Finally, we use the definition of quantum mutual infor¬ 
mation 


S(Pdmm) ~ ${pd) + S(pmm) I(D : MM) 

and the fact that A S B = 0, again from the periodic 
steady state condition on D, in Eq. (B3) to obtain Eq. (6) 


Qh^c(Pc fih) + MM — ^ P D:MM 


D(Pdmm\\(^dm®Pm) D(p D mm\\p^dm®Pm^ — 0 (Bl) 
where 


Equation (4) follows in the same manner, but without 
including M, and Eq. (3) requires neglecting correlations 
(A Id.m ~ 0). 


D{p\\v) = S(p) ~ tr(plnd) 


is the quantum relative entropy. Using 


tr 


In (p A <S> Pb ) = In pa 0 t B + t A ® hip# 
Pab ( Oa ® = tr ( paOa^ = ( Oa ) 


Appendix C: Matrix Product Density Operator 
Update Method 

We describe our matrix product quantum operation 
shown in Fig. 2(a) by its action on the product basis 
elements in Eq. (7). Using a generalization of the result 
in Ref. [58] from vectors to operators as in Ref. [49], we 
update each element by n sequential interactions, 4> nr , 
as 


and 


tr 


(Pd ~ Pd)^Pd P) 


= 0, 



y]tr { B lD V q C l P 


3 = 1 


n 



3 = 1 


where 
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(&’*) 

| = tr 

bn ® VmY <t>T Gn ® °m) 




- 

(«“■) 

| = tr 

a/3 

(°d)Ad{ o) 



where 5i D p is the Kronecker delta and the a lD are the 
basis elements for D. These expressions follow from com¬ 
pleteness under the trace inner product for this operator 
basis, 


After the next interaction, we have 


Pdmm 


N 

tr [ <g> i) P JJ 




j=n +2 



(C3) 


N 



j=n +1 


A,B) = trA^B 

/ tr 


Thus, after rearranging, 


N 

pDM{nr) = y^tr ^B %D 0 Aj V (t 


We see from this expression that its convergence with n 
depends on that of the transfer matrix (^2 k D k 0 A k ^j to 

its fixed point. It is therefore exponentially convergent at 
the rate of the inverse correlation length of the outgoing 
MPDO, as we expect. 

We now simplify Eq. (C2) using is the classicality of 
our interaction, </> T , which implies = 0. Rearranging, 
and making the identifications 


x^ni 

( V C^ k ® A k ) 

N 

A £, -5> 

| 

fn**') 

j=i 

\ fc / 

j=i 



Vi- 1 /J 


a 


id 
D > 


(Cl) 


and in the next interaction, we make the replacement 


~(k) 

PkP M M 




* (n A 

\j=n+1 ) 



N 



j=n +1 


i 0 A in+1 ^2 c in+lk 0 A k . 

k 


We see that this leaves the description of Eq. (Cl) ap¬ 
proximately unchanged for sufficiently large AT, n. This 
is the aforementioned periodic steady state of the system. 

Finally, we trace out all previously interacted degrees 
of freedom to obtain 


Pdmm 


N 


id^ 


j=n +1 



N 



j=n +1 


(C2) 


for k G {0, l} n gives our Theorem 2, 


where 


E q/c) _ q/c) 
PkPD ® PMW 

k 


p§ 

j=i 


Ak) 

0 ~ 
MM 



M 



Pm(0) 


and 


where 

D k = ^2C lk tia l . 

i 


T kj ( PD ) = tr M 4>r (pD ® <r^) , 
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The classicality of the p^ follows from the uniqueness 
of the fixed point of the transfer operator T k , which is the 
periodic steady state of the demon upon interacting with 
a string of uncorrelated pure classical qubits, \k)® n , and 


found in Ref. [ ]. Because this fixed point is classical and 
unique, T k must be dephasing in the energy eigenbasis, 

and so the are classical as well. 
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